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Introduction 


It has been demonstrated in [6], [7] that the Legendre-tau approximation 
is a quite powerful approximation method for hereditary differential systems 
in many instances. However, there remained an important question which had 
not been resolved. This is the question concerned with the presentation of 
exponential stability under approximation. We establish it in Section 5. 

As observed in [7], the Legendre-tau approximation scheme provides a good 
approximation technique for computation of eigenvalues for hereditary 
differential system as well as its optimal closed-loop system. We give a 
justification of this observation deriving a characteristic equation for 
eigenvalues of the approximating system and relating it to the Pade 
approximations of the exponential function. Moreover, it leads to a 
characterization of detectability and stabilizability conditions for the 
approximating system and the preservation of those properties under 
approximation. 

The results discussed in this paper are similar to those for the 
"averaging" approximation scheme [1] that have been obtained in [1], [5] and 
[9], and a great deal of our discussions are motivated and inspired by those 
investigations. We refer to [1], [8], [2], and [5] for the the summary of the 
earlier contributions on optimal control and numerical approximation problems 
for hereditary differential systems. 

The following is a brief summary of the contents of this paper. In 

Section 2 we state the type of problems to be considered and review the 
equivalence results between hereditary differential equations and abstract 
Cauchy problems on the product space I? x and then describe the Legendre- 
tau approximation scheme within the abstract framework. The convergence of 
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the ad joints semigroups S N (t)* is proved in Theorem 2.2. In Section 3 we 
derive a characteristic equation for the approximating system and discuss its 
stability and convergence properties. In Section 4 we review results on the 
linear quadratic optimal control problem and prove the convergence of the 
solutions of the approximating evolutional Ricc.ati equation in trace norm. In 
Section 5 we prove the Gibson conjecture; i.e., that exponential stability is 
preserved under approximation. An important consequence of this is that by 
following the approach given in [5] one can establish the convergence of 
solutions of the approximating algebraic Riccati equation in trace norm (which 
implies strong and uniform convergence). 

The notation used in this paper is standard and exactly the same as used 
in [6], [7]. We denote by Z the product space I? x I^Cl-r.O]; I^ 1 ). Given 
an element z e Z, q e I? and $ e denote the two coordinates of 

z:z = (q,$). For any function $ of the independent variable 0, we shall 

* 9 

use $ or <|> to denote the derivative of <f> with respect to 8. 

2. Legendre— Tau Approximation 

In this paper we will restrict our analysis to the system: 

d 0 

x(t) = Agx(t) + Ajx(t-r) + / A(0)x(t+0)d9 + Bu(t) 

" r 

x(0) = q, x(0) = <f>( 0 )» -r _< 0 < 0, 

where 0 < r < +» and Aq, Aj and A(») are nxn matrices, the elements of the 

latter being square integrable on [-r,0]. It is well known [1] that for 
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,N-1 3 N 


3 N, 


t z (t ,9) = “ z “(t,8) 


d N 


z“(t,0) = A Q z l, (t,0) + A lZ w (t,-r) + / A N (9)z N (t,0)d6 + Bu(t), 


-r 


N N-l 

where A (0) = (P A)(0), - r < 0 < 0. As shown in [7], using (2.3), 
z N (t), t > 0 satisfies 


where 


3 N, . 

u z (t) 


,NN. , „ , . 

A z (t) + Bu(t) 


z N (0) 


Q N (n,<J>) 


. N N N « T v i 

A = QAL, N>1. 


(2.4) 


N a' t 

If we denote S (t) = e A , t 2 0, then it is proved in [7] by using the 
Trotter-Kato semigroup approximation theorem [10], that S N (t)z converges 
strongly to S(t)z for t ^ 0 and z s Z, and the convergence is uniform 
in t on bounded intervals. 


Lemma 2.1. For (y,i|>) e Z, 

<A N ) (y.'lO = (A^y + Ao), - + A N («) T y) e z, 


N N-l , T 
* - P + (Ajy - (P 


N-l , 


'l')(-r))p N . 


where 
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Proof : Note that for e Z, 

A N (n ,<|>) = (A Q n + A^ N (-r) + J° A N (9)<(.(0)d0 , * N ) 


where 

^ = pN-1^ + (t, _ (P N_1 (),)(0)) p N . 


(2.5) 


Then we have 

' I = <A N (n,«j,),(y,t)> z . 

= <A Q n + A^Vr) + / A N (0)<K0)d0,y> 

-r 

+ I <f (0) ,^(0)>d0 

-r 

“ <n »Agy> + / <<f»(9), A N (0) T y>d0 

-r 

+ <* N (-r), A*y> + <n,/(0)> - </(-r),/(-r)> 

- J° <<|> N (0) . ^ N (0)>d0, 

-r 

where $^(0) = n and is a vector in if 1 whose elements are polynomials 

of degree less than N-l. Note that = A^y and P N *<|> N = $ N . Thus, 

I = <n,Agy + i}» N (-r)> + / <$(0),A N (0) T y - ^ N (0)>d0, 

-r 

for all (n ,<!>)> (y ,4») e z, which completes the proof. (Q.E.D) 
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Although we will not pursue the details here, the error estimates in [6], 
[7] can be used to prove that 


M * 

(A ; (y,*) 


A (y,t), 


for all (y,<f>) e £>f(A ) 2 ), 


Hence the Trotter— Kato theorem again applies to obtain 


Theorem 2.2: 


N * 
S\t) z 


S (t)z. 


t > 0 


uniformly on bounded t-intervals. 


3. Characteristic Equations 

In this section we will derive a characteristic equation for If 

N n ' 

z e Z is an eigenfunction of A • 


Nx N 


(Al - A )z = 0 for some complex number A, 


(3.1) 


N 

then z = 


n N N 

Q z 


N-l 


N A N 

(h > 1 a i. Pu j ** L et L^z = (n N ,<f> N ) e V(A) where 
k=0 


N-l 

.N N N 

♦ -J 0 V k + V„ 


N-l 


N N r N 

N - n - I a k * 


k=0 


Then from (3.1) and the definition (2.5) of A N » we have 
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From (3 


where 

by 


for N 
define 


Xn N - A Q n N - A 1 <|) N (-r) - / A N (0)<J> N (0)d0 = 0 

-r 

x(p n -V) - * N = 0. 


r n I 


r N “I 

a o 


a o 

• 

2 „ 

• 

• 

- — S 

• 

• 

r 

• 

N 


N 

_ a N-l_ 


L a Nj 


S is the matrix representation of the derivative (3/30) and 



0 

1 

0 

1 

0 

• • 

1 

0 


0 

0 

3 

0 

3 

• • 

0 

3 


0 

0 

0 

5 

0 

• • 

5 

0 

S = 

• 

• 

• 

• 

• 





• 

• 

• 

• 

• 

• 

• 

• 

• 

• 





0 

0 

0 

0 

0 

• • 

2N-3 

0 


_ 0 

0 

0 

0 

0 

• • 

0 

2N-1 _ 










even. For 

N odd. 

only 

the last column 

of S 

is different 

r K N, N 

{ V k=0 

by 










b? - 

h? , 

N 

= a, > 

0 1 k 1 

N-2 





k 

k+2 

k’ 







, N 

N 


l n n 






b , 

= a„ 

j and 

b.f B ci.. • 






N-l 

N- 


N N 





(3.2) 

(3.3) 

(3.4) 
is given 

Let us 


Then from (3.4) if 


W - (f)A 


where 


Ub 


1 


. N 

0 


b l 


m 

• 

m 

• 

= e ( y ) 

• 

• 

• 

0 


b ; 

— 


_ N_ 


E 


N 


is a tridiagonal matrix: 


e n (p) 



y 

3 

-y 



Applying Cramer's rule, one obtains 


where 


,N k N ,N, N 

\ ■ “ Vi b o /lf i 


11 , 
Y k - det 


2k-l 

-P 


-y 2N-1 J 


, 0 < k < N and y 


N+l 


Since Y^ = (2k-l)Y^ + ^ + y^ Y^> it follows from (3.5) that 


N / 01 v k N N . N 
\ = (2k+l)p T k+1 b 0 / Tl , 


0 < k < N-l 


N N ,N, N 

a N =U V Y 1* 
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M K7 VT 

Note that n - £ af, - bj! + b.. It now follows from (3.2) that 

k=0 k U 1 


(XI - A 0 )(yN + yy 2 ) ~ “ MY 2 ) “ / A N (0)5 N (0)d8 = 0, 

-r 


where 


5 N (9) 


N 

- I 

k=0 


k N 


(2k+l)y Y 


k+2 ^k 


/’28+r-v, N 

p u( r J b Q 


and 


Thus, we obtain 


.v - (£)*• 


(u- 


N N 0 
Yi > UY 2 


A — A 

0 1 N ^ N 

Yi + bY 2 


- / A N (0)<|> N (0)d0)bg = 0 (3.6) 


-r 


where 


* N (0) = l (2k+l)u k Y l +2 / (Yj + MY 2 ) P k (^ £ ), -r < 0 < 0. 


k=0 


For N > 1, let us define 


.N N , 0 N N 

6 = Y^ and 0 = Y 2 « 


Then we have the recurrence formula: 
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6 N = (2N— 1 ) 6 N_ 1 


+ U 6 


2 -N-2 


(3.7) 


N N-l 

3 = (2N-1)8 


. 2 .N-2 

+ u 8 


with 6^ =6^ = 1, 3^ = 0 and 8* = 1, so that if = 6^ ± y8^> then the 
p's satisfy the recurrence formula: 


N 

P 


(2N-1 )p 


N-l 


+ V 


N-2 

P 


with p^ - 1 and p^ = 1 ± y. The recurrence formula (3.7) is exactly the 
same as that for the diagonal Pade approximation of exp(-rX) [4] . Note that 


N 

<|> N (0) = 1 and <J> N (— r) = — ^ = Pade (e” rX ). (3.8) 

P+ 


Hence from (3.6), we obtain: 


Theorem 3.1 X is an eigenvalue of A N if and only if X satisfies 
det A^(X) = 0 where 

A N (X) = XI - Aq - A, Pade (e" rX ) - J° A N (0)<f> N (0)d0. 

-r 


Lemma 3.2: If the elements of A(») are absolutely continuous on [-r,0], 

then there exists a positive constant C such that 

|/ A N (0)<(. N (0)d0| < C 

-r 
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for every X e £ with Re 1 > 0. 

N 

Proof ; Since <j> . satisfies 

N • N N 

H - = XaV)p N 

N N N 

where a (X) = y /p + » we have 

,N, 0 . X9 . T °. N,,. X(9-0 r 2C+r bc 
$ (0) = e + J X a (X)e ’'p Jd?. 

0 

Thus , 

J A N (e )<*> N (e )d0 - J A(0)<j> N (0)d0 - J A(9)a N (X)p N d0 
-r -r -r 1 

= J° A(0)|e X9 + /° Xa N (X)e X(0 “ 5) p N d5 - a N (X)p N (^^)j d0. (3.9) 

Here, 

1° A<9> I 0 x e x < e -«)p N (Bii)d 5 
-r 0 

- /° i A(9)Xe i(e - 5) de P N (^)d5 

-r -r 

- J° /a( 0 - A(-r)e _X(r+ ° - J 5 e X(9 “ C) A(0)do\ P N dC. 


Hence from (3.9) 
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0 0 0/5 
1/ A N (0H N (0)d9| < / |A(6)|d0 + |a N (X)| / (|A(-r)| + / |A(0)|d0 

-r -r -r \ -r 

0 / 0 

</ | A(0 ) | d0 + r |A(-r)| + J | A(0) |d9 

-r V -r 

on (J: + = {X e <£, Re X >. 0} where we used |a N (X)| _< 1 on <|: + . (Q.E.D.) 

Corollary 3.3: If the elements of A(*) are absolutely continuous on [-r,0] , 

N 

then the complex function det A (X) cannot have a zero In the closed right 
halfplane outside the disc of radius |Aq| + | Aj | + C. 

Proof ’ Note that |Pade (e rX )| _< 1 on <t + . The corollary follows from 
Lemma 3.2. (Q.E.D.) 

Lemma 3.4: A N (X) converges to A(X) uniformly on every bounded subset of 

the complex plane . 

Proof : We only need to show that 

0 „ „ 0 

/ A (0 )<{> (0)d0 ► / A(9)e A °d9, uniformly. 

-r -r 

•From (3.9) 

0 0 

|/ A N (0)<|> N (0)d0 - / A(0)e A8 d0| 

-r -r 

= |/ a N (X)A(0)^/ Xe X(0 " O p N d5 - p N (^il)jd0| 
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|/ a N (X) 


-r 


-A(5) + / Xe X(0 " 5) A(9)d0 


-r 




< |a B (X)|/l + r/° |Xe X6 | 2 dsj 1/2 |A(6)| 2 dej 1/2 ' /2 , 


which converges uniformly to zero on every bounded subset of the complex 
plane. (Q.E.D) 


4. Riccati Equations 

Let G be a non-negative, self-adjoint operator on Z and C be 

a pxn matrix. Consider the optimal control problem on a finite interval: for 

given initial data z = (n , <J> ) e Z 

minimize J(u; [0,T]) = / (|Cz(t)| 2 + |u(t)| 2 )dt + <Gz(T),z(T)> (4.1) 

0 

over u e LgUO,!]; 1?“) subject to (2.2), where C (n »4> ) = Cn for 

(n,<|>) e Z. It is well known [3], [5] that the optimal solution u° to 

(4.1) is given by 

f) * 0 

u U (t) - - B n(t)z (t) 

where II (t), t < T is the unique non-negative, self-adjoint solution to the 
Riccati equation: 

4 - <H(t)z,z> = -2<Az,n(t)z> + <B*n(t)z,8*n(t)z> - <Cz,Cz> (4.2) 


for all z e V(A) 
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tl(T) = G, 


and z^(t) satisfies the evolution equati 


on 


z°(t) = (A- B8*n(t))z°(t) 


z (0) = z. 


Consider the Nth approximate problem to (4.1): minimize 


J N (u;[0,T]) = / ( |C z N (t) | 2 + |u(t)| 2 )dt + <G N z N (T),z N (T)> (A. 3) 

0 


subject to (2.4) where G^ = Q^GQ^. The optimal control to (4.3) is 

given by 

N, \ * _N, N N, . 

u (t) = -B n (t)z (t), 

where II^(t), t < T is the unique non-negative, self-adjoint operator to the 
Nth approximate Riccati equation: 


dN NN *N *N 

^ <n (t)z,z> = -2 <a z,n ( t )z> + <b n (t)z,B n (t)z> - <cz,cz> 

for all z e Z (4.4) 


At) = g n 


and Z N(t) satisfies 
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d_ 

dt 


z N (t) = (A N - 88*n N (t))z N (t), 


t > 0 


N 

z (0) = Q 


N 

z. 


From Theorem 2.2 and Theorem 6.1 - 6.3 in [5], we have the desired convergence 
results: 

Theorem 4.1. n N (t) converges strongly to II(t) for t _< T and the 
convergence is uniform for t in bounded intervals. If u^ = -8 ll^(t)z^(t), 

A VT 

t < T where z w (t) is the mild solution to 


^?(t) = (A - B8*n N (t))z N (t) 
z N (0) = z; 

i.e., the Nth feedback -ontrol law applied to the original hereditary system 
(2.1), then for any e > 0 there exists a nondecreasing function 
N_ (• ) : [0,®) *■ IR?' such that for n > N_(T) 

j(u N ; [0 ,T] ) _< J(u°; [0.T]) +elzll 2 . 

Moreover, if G is given by G(n,<|>) = (G 0 n,0), then n N (t) converges in 

trace norm to H(t) for t < T and the convergence is uniform for t in 
bounded intervals. 

Let us now consider the optimal control problem on the Infinite interval. 
For given initial data z = (h,<f>) e Z, minimize the cost functional: 
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OO 

J(u,z) = / ( | Cz(t) | 2 + |u(t) | 2 )dt 

0 


subject to (2.2). 


(4.5) 


Definition : The pair (A»B) is stabilizable if there exists an operator 

K e Z. ( Z , I?) such that A — BK generates a uniformly exponentially stable 
semigroup. (C,A) is detectable if (A*,C*) is stabilizable. 


The following theorem is now standard [11], [5]. 


Theorem 4.2. _If (A , 8) is stabilizable and (C,A) is detectable, then the 
algebraic Riccati equation (ARE) : 

(All + IIA — IIB8 II + C C)z = 0 for all z e 17(A) 

has a unique self-adjoint, non-negative solution. Moreover, if II denotes 
the said solution, then A - 88 II generates a uniformly exponentially stable 
semigroup and the optimal control to (4.5) is given by 

0 * 0 

u U (t) = -B IIz (t) 

where z^(t) is the mild solution to 

ij_ z °( t ) = (A - B8*n)z°(t) 


z^(0) = z 
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Remark: 

(i) 


(ii) 


(iii) 


The 


For the hereditary differential systems, we have 


(A,B) is stabilizable if and only if 


rank[A(X),B] = n for all X e <jl + , 


(C,A) is detectable if and only if 

= n for all X c 
a complex number X e o(A — 8B H) if and only if Re X < 0 and X 

a 

satisfies det A(X) = 0 where 




rank 


A (X ) 


A (X ) = XI 



-BB 


T I 


-D(-X) T 


D(X ) = A Q + e" rX + / A(9)e X9 d0[5]. 

-r 

next lemma gives the characterization of stabilizing and detectability 


for the approximating system (2.4) which is precisely the analogon to those 
given in Remark. 



-19- 


Lemma 4.3. 

( i ) ( > B ) is stabilizable tf and only if 


rank[A N (A),B] = n for all A e <J; + 


(ii) ( C, A n ) is detectable if and only if 


rank 


A N (A) 


= n for all A e <£ + . 


Proof_: We will only prove the statement (ii) since (i) follows from Lemma 2.1 
and (ii) by duality. It follows from the finite dimensional linear system 
theory that (C,A N ) is detectable if and only if 


A^7. = Az and C z = 0 for A e 

imply Q^z = 0. Hence from (3.6), (3.8) and Theorem 3.1, it is equivalent to 

ker A N (A)Oker C = {0} for A e <J: + . (Q.E.D. ) 

The next corollary follows from Corollary 3.3, Lemma 3.4 and Lemma 4.3. 

Corollary 4.4. Suppose that the elements of A(*) are absolutely, continuous 
on [ —r >0]» If (A,B) is stabilizable ((C,A) is detectable, respectively), 
then (A N , B) is stabilizable ((C,A^) is detectable, respectively) for N 
sufficiently large. 
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For the rest of this section we assume that the conditions stated in 
Corollary 4.4 hold. Let us consider the Nth approximate problem to (4.5): 

minimize 

J N (u,z) =/ (|Cz N (t)| 2 + |u(t) | 2 )dt (4.6) 

0 

subject to (2.4). The optimal control u N to (4.6) is given by 


N, . *N N, , 

u (t) = - B n z (t) 


and 


N * N 

N, * (A - 88 n )t n N 
z (t) = e v Q z, 


where IT** is the unique nonnegative, self-adjoint solution to (ARE)* 


(AV n N 


+ n N A N - n N bb* 


n N + c*c = ({). 


In the following lemma we give a characteristic equation for eigenvalues of 

,N * _N 

A - BB n . 


Lemma 4.5: 

A complex number 

X is 

an eigenvalue of A^ 

if Re X < 

0 and X satisfies det 

A N (X) = 0 where 



d n (X) 

T 

-BB 


I 

M 

II 

*< 
< < 





T 

N, . ' 



-C C 

-D (-X) 


N _ ttN 


D N (X) = A 0 + A 1 <f» N (— r ) + / A N (0)<f> N (e)de, 
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Proof. It follows from [5] that \ e a (A N - 88* II N ) if and only if 
Re X < 0 and X e a(H^) where 


-88 


h n = 


* 

■cc 


N * 
(A ) 


on Z N xZ N 


and Z N = Q N Z. 

From (2.5) and Lemma 2.1, if X e a(H^), then there exists an element 
((h,?), (y,t)) e Z N x z N such that 

A Q n + A^(-r) + / A (9)<K9)de - BB y = Xn 
r 

l = X$ 

-C T Cy - Aq y - tp(0) = Xy 

NT • 

A (•) y - if> = X\p 


where 

<|> = 4> + (n - <|)(0))p N , n = <f>(0) 

'I' = 'I' + (A 1 y - T|;(-r))p N and ij>(-r) = y. 


Thus the similar arguments as given in Section 3 allow us to conclude that 
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X e o(H N ) if and only if det A N (X) = 0. 


(Q.E.D. ) 


5. Gibson's Conjecture 

In this section, we prove the Conjecture 7.1 in [5] for the Legendre-tau 
approximation: if the semigroup {S(t), t > 0} is uniformly, exponentially 

stable, then there exist positive constants M and w such that, for N 
sufficiently large 

ns N (t)a < . 


To this end we need the following two results . 

lemma 5.1: Let us denote by Aq, a generator on Z defined by P(Aq) = P(A) 

and 

A 0 (4>(0),<j>) = (-*<(>),♦) e Z 


and 


define A^ = Q N A q L N , N > 1. Then there exist positive constants ^ 


and Mo such that 


A o t "V 

Je I < M Q e 


and 


* ^ '2 , M1 ,2 

jf M o" z " 


/ | E e z 

0 


where E is an operator on Z defined by 



E(n,<t») = <}>(-r) 


for (n,i}>) e Z. 


Proof: We define the weighted inner product <•.•> on Z bv 

g 


0 


<(n ,*),(y ,*)>„» <n,y> +/ <*(0M(6)> g(9)d9 

R n 


g " R n -n 


(5.1) 


where g is positive on [-r,0] defined by 


3 j. 1 r 20 + 
g< 0 ) = 7T + 7T ( — • >' 


Since -j <!. g(9) _< 1 on [-r,0], if Zg denotes the completion of Z with 

respect to the inner product (5.1) and 1*11 denotes the induced norm, then 

g 

IzI 2 < tz« 2 < 2!zl 2 . (5.2) 

g - - g 


For z N = (n,4») G Z N = Q N Z 


T . N N N. r ._IN-1 ,JN TIN. , 0 . I i. 

I = <A 0 Z ,z > = J <P <J> >g(9) - |n| 

s -r 


0 


„N-1 . N *N V 


where <j> N » $ + (n - $(0 ))p n and <|> N (0) = n. 
Note that 


/ <n-<K0))P N (— * N > g(9)d9 


-r 


= / <(h - <K0))P N , ft £ (^-t£)d9 


-r 
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where we 
NP n (x) = 


I 


where u) 


= / |n - K0)| 2 p n p n i 

-r 

1 1 ? 

.= 4 "/ |n - K0)| P N (x)P N (x)x dx 

1 1 9 

= i-/ (2N— 1) |n - <KO>r P N (x)P N _ 1 (x)dx 

-r 

1 1 2 

" 4 J In - *(0)1 P N (x)(NP N (x) - (N-l)P N _ 2 (x))dx 
-r(wr) ln-W,| 2 , 

N— 1 • M • M 

used that P <{>=<)> and 

(2N-l)xP^(x) - (N-l)P^_ 2 (x) » “1 <_ x 1. Thus, we have 

- / <* N ,i N >g(e)de - |n| 2 - j |n - M0)| 2 

-r 

_ 1 |,N/ n > | 2 1 ,.N, x | 2 1 ,2* f ° | N | 2 JO 

B 2 !♦ (0) l " T 1+ (-r) l ~ s' ( 7> / I* I d9 

-r 

" I^l 2 -i C-^rJ ln-<K0)| 2 

- -I H 2 -J lA-r)| 2 -i (1) J° M 2 d9 -i In - *( 0 )| 2 

-r 

< - ai 0 i(n ,<!>)» 2 | (— r ) | 2 , 

j = min^ , . It now follows that 

N. . N, .. N. . 2 1 N>. v|2 

Z Z ( t )> g £ “ “o” Z ” IT I E Z 
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or 

Bz^(t)B 2 < Bz^(0)B 2 - 2 oj / Jz N (s)H 2 - -i- / |Ez^(s)| 2 ds (5.3) 

g- g 0 q g 4 q 

< B z N (0)n 2 - 2u)„ / !z N (s) 1I 2 ds. 

- g 0 o g 

Hence by Gronwall's lemma 

N -a> 0 t N 

Bz (t)B < e Iz (0)1 , t > 0. 
g - g - 

From (5.3), for all t _> 0 

J t |Ez N (s)| 2 ds < 4Bz N (0)B 2 < 4Bz N B 2 < 4BzB 2 . (Q.E.D. ) 

0 g 

Corollary 5.2. Let I denote an operator on Z defined by I(n,<|>) = n e S 1 
for (n,<|>) e Z. Then, If 

00 N 

j \le t z| 2 dt £ KBzS 2 
0 

for some positive constant K, there exists positive constants coi and Mj, 
such that 

,N -w t 

«e At l <Mj e , t > 0. 

Proof: The same argument as in the proof of Lemma 5.1 yields that for 

(n,+ ) e z N 
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<A^z^,z^> = <AqT 1 + Aj<J>^(-r) + / (9)4>^(9)d6 ,n> 

g -r 

+ 7 hi 2 lA-r)! 2 -£<i> 1 ° |K0)| 2 d8 

-r 

<(£ + |A 0 1 + |A *J| +2r/° |A(6)| 2 d0) |n | 2 - j u) Q IzV 

-r 

II 2 1 , N B 2 

= a |n | - j t * 0 tz n , 


where we used the relation: 2 <x,y> j< |x| + |y| for x,y e IP. It then 

follows that 

< dt Z Z ( fc )> g 1 a H ) I “ 2 “q " Z 


or 


n z N (t )n 


2 

g 


+ UK 


/ Jz N (s)H 2 ds < n Z N (0) B 2 + 2 a / 
0 " g 0 


| n N (s ) J 2 ds 


for all t > 0. 


Hence we obtain 


/ Hz N (s)R 2 ds <^ — (1 + 2o)Rzl! 2 . 

0 “o 


The corollary now follows from Datko's theorem (see Lemma 7.4 in [5]). 


(Q.E.D.) 


We are now ready to state the main theorem. 
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Theorem 5.3: If the elements of A(») are absolutely continuous on [-r,0] 

and the semigroup {S(t), t ^ 0} Is uniformly, exponentially stable, then 
there exist positive constants w and M such that, for N sufficiently 
large 


,N 

, A t. . „ -tot 
Be B < Me 


Proof : It follows from Corollary 3.3 and Lemma 3.4 that there exists an 

Integer Ng such that, for N >_ det A^(A) * 0 for all A e $ + . But for 
M >_ |Aq| + |aJ + C = 0 and u e 1, it follows from Corollary 3.3 that 


|A N (im)" 1 | < 1 


- m - B 


This inequality, when combined with Lemma 3.4 shows that there exists a 
constant a such that for N > 


A N (iaj) _1 1 <o, Bel. 


Note that for z e Z 


t A 0 t A N (t-s) 

z = e z + J e I 


f (s)ds 


N. 


(5.4) 


where f(t) = Fe z, t > 0 and F is an operator on Z defined by 


F(n,«J>) = (I + A 0 )n + A^c-r) + / A(0 )<j>(8 )d9 for (n,4>) e Z. 


-r 
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It then follows from Lemma 5.1 that there exists a positive constant Y such 
th8t 

°o a t 00 _ 

/ | X e ^ z|^ dt and / |f(t)|^ dt ^ yBzH . 

- 0 0 

t I N/. \ ^ A 

If y(t) = J 1 e A K ~ J I f(s)ds, then the Fourier transformation y(iw) of 
0 

y is given by 


y(iw) = A N (i<o) * f(ito). 


Hency, by the Parseval's equality. 


/ | y ( t ) | dt = / | y (ico) | d(o 

0 


w w 

<a 2 J | f (ico) | 2 dco = a 2 J |f (t ) | 2 dt. 

—OO 0 


It now follows from (5.4) that 


00 N 

( I Ie A 1 z 1 2 dt < 2y(l + a 2 ) 0 zE 2 


0 


which completes the proof along with Corollary 5.2. 


(Q.E.D. ) 


The following theorem is an important consequence of Theorem 5.3. 

Theorem 5.4. If the elements of A(») are absolutely continuous on [— r , 0] , 
rank C = n and (A, 8) is stabilizable, then for N sufficiently large 
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Remark : In the above proof the condition: rank C = n is only used for the 
derivation of (5.5). But it seems to be enough to assume the detectability of 
(C, A) for such a derivation. 



-32- 


References 

[1] H. T. Banks and J. A. Burns, "Hereditary control problems: Numerical 

methods based on averaging approximation," SIAM J. Control 
Optimization, 16 (1978), pp. 169-208. 

[2] H. T. Banks and F. Kappel, "Spline approximations for functional 
differential equations," J. Differential Eq., 34 (1979), pp. 496-522. 

[3] R. F. Curtain and A. J. Pritchard, "The infinite dimensional Riccati 
equation for systems defined by evolution operators," SIAM J. Control 
Optimization, 14 (1976), pp. 951-983. 

[4] B. L. Ehle, "A - stable methods and Pade approximations to the 
exponential," SIAM J. Math. Anal., 4 (1973), pp. 671-680. 

[5] J. S. Gibson, "Linear-quadratic optimal control of hereditary 
differential systems: Infinite dimensional Riccati equations and 

numerical approximations," SIAM J. Control Optimization, 21 (1983), pp. 
95-139. 

[6] K. Ito and R. Teglas, "Legendre-tau approximation for functional 
differential equations," NASA CR-172155, ICASE Report 83-17, NASA 
Langley Research Center, Hampton, VA, June 1983. 



-33- 


[7] K. Ito and R. Teglas, "Legendre-tau approximation for functional 
differential equations, Part II: The linear quadratic optimal control 

problem," NASA CR-172397, ICASE Report 84-31, NASA Langley Research 
Center, Hampton, VA, July 1984. 

[8] A. Manitius, "Optimal control of hereditary systems" (lecture notes), 
in Control Theory and Topics in Functional Analysis, Vol. Ill, pp. 43- 
128, International Atomic Energy Agency, Vienna (1976). 

[9] D. Salamon, "Structure and stability of finite dimensional 
approximations for functional differential equations," MRC Report 
#2586, University of Wisconsin, Madison, WI, October 1983. 

[10] A. Pazy, Semigroup of Linear Operators and Applications to Partial 
Differential Equations , Springer-Verlag, 44 , 1983. 

[11] R. B. Vinter, "Filter stability for stochastic evaluation equations," 
SIAM J. Control Optimization, 15 (1977), pp. 465-485. 



1. Report No. NASA CR- 172447 
ICASE Report No. 84-42 


2. Government Accession No. 


3. Recipient’s Catalog No. 


4. Title and Subtitle 

LEGENDRE-TAU APPROXIMATION FOR FUNCTIONAL DIFFERENTIAL 
EQUATIONS, PART III: EIGENVALUE APPROXIMATIONS AND 

UNIFORM STABILITY 


5. Report Date 

August 1984 


6. Performing Organization Code 


7. Author(s) 

Kazufumi Ito 


8. Performing Organization Report No. 

84-42 


10. Work Unit No. 


9. Performing Organization Name and Address 

Institute for Computer Applications in Science 
and Engineering 

Mail Stop 132C, NASA Langley Research Center 
Hampton, VA 23665 


11. Contract or Grant No. 

NAS1-17070 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 


13. Type of Report and Period Covered 
Contractor Report 


14. Sponsoring Agency Code 

505-31-83-01 


15. Supplementary Notes 


Langley Technical Monitor: Jerry C. South, Jr. 

Final Report 


16. Abstract 

The stability and convergence properties of the Legendre-tau approximation for 
hereditary differential systems are analyzed. We derive a characteristic equation 
for the eigenvalues of the resulting approximate system. As a result of this 
derivation we are able to establish that the uniform exponential stability of the 
solution semigroup is preserved under approximation. It is the key to obtaining the 
convergence of approximate solutions of the algebraic Riccati equation in trace norm. 


17. Key Words (Suggested by Author(s)) 

hereditary difference systems 
Legendre-Tau approximation 
Pade approximation 


19. Security Classif. (of this report) 
Unclassified 


18. Distribution Statement 

64 — Numerical Analysis 
Unclassified - Unlimited 


20. Security Classif. (of this (sage) 
Unclassified 


21. No. of Pages 
35 


22. Price 
A0 3 


N-30S 


For sale by the National Technical Information Service. Springfield. Virginia 22161 













